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/ I. INTRODUCTIOrf AND SUMMARY ' ♦ 

The purpose of this study iX, to' invesii^ate "the use of the Defense 
Mapping. Agency Aerospace Center (OMAAC) <mital Source .data for simulation 
•of Forward Looking Infrared (FLIR) and Low' Light Level Television' (LLLTV) 
sensor» imagery during various periods of-thi day and s'easp^ of the year. 
This includes determination of those deficiencies of the DMAAC data ba.se 
which -limit, its use as well as 'identification of the additional parameters 
which would -increase its utility for sensor simulation. In addition, the 
.effort involves development' of Computer Image Generation (CIG) algorithms, 
to simulate the effects of diurnal and seasonal variations on sensor imagery 

In.order to fulfill the intended goals of this study, the approach 
_jrfas. divided into four tasks which are described in the following four" 
sections gf this report. The first task involved construction of a data 
base for a chosen target ar^a for three different resolutions of the DMAAC 
data base. The subject target area was one^of 16 available from the AFHRL 
Seasonal Sensor Handbook: The second task consisted of the- derivation of < 
fairly sophisticated visual a;id" thermal tonal models from thfebretical 
thermo-electromagnetic models. These models blend the physical properties 
of the scene withe the- environmental conditions and particular sensor char- 
acteristics to result in an accurate simulation of sensor imagery under 
various diii<;nal and seasonal conditions. The third task consisted of 
generating siWlated sensor imagery for various diurnal and seasonal .|| 
conditions using the data bases and tonal models developed in the first 
two tasks. The fourth task consisted of. analysis of the deficiencies af 
the simulated imagery, in order to determine what parameters need to be 
added to the DMAAC data base to support 'effectWe Simulation of sensor 
imagery. 



At the beginning of the program the intent was .to. generate the imagery 
witli a modest extension of the ^ata, base requirements. However, some Short- 
comings were assj^ciated with thi"S approach* The 12 material cades used in' 
* the DMAAC- specif ication do not fully^-rpflect the physical progertjes of the 
^material- For ex-ample, o'n' the sceri^evaluated inrthis study, both the 
cement p^lamt/ -and the hospital -illve Uhe 'same' materia.1 cocle'yet they have 
difference ihennSl propfeuties.- A^so, the DMAAC feature code is use oriented 
and is npt related to the reflective or emissive properties of the buildings 
Thus, the hospital has a feature code that is identical to- that of schools' 
and' churcnes, all of which are identified as institutional structures, -but' 
this identification has nothing to-do wi^h the''^underlying physical 
'characteristics. With limited additional descriptors j'fi t^e data base, it 
was felt that only a little more fidelity could' be added to the simulated 
passive imagery which is derived throiigh transformations of the DMAAC data. 

A basis of much of t;he tonaVmod^ling is theVesul^tJttj^^^ 
program under Air Force Human Resources, Laboratory (AmBw^^ * 
which resulted in'the AFHRL Seasonal Sensor Handbook. "^j^jfj^jB^^ments 
were ta-ken on multiple days in each'of the four s^^s^ns Q^^mB^lk^r^ on 
16 different target scenes. yThe instrument and ehvironmental^ohditions 
were recorded. THis handbook provided a qual itativ^ /source during the 
formatien of the tonal models. ^ Also, once , the models and i'mageryt^Ve 
derived, the handbook v/as u^sed for qualitative evaluation. ^ - 

Thus, the approach taken in the study was to determine what information 
is needed to simulate the imagery contained in the handbook. Then imagery 
.using lesser amounts of data" can be generated^ for various user applications. 
The impact on the data base can then be assessed for a given application 
and level of detafl in the 'imagery. * . - 

.-./■• ^ . 

The results of tMs study i/iclude a realistic tonal model for the 
simulation of passive sensor^.imagefy of ^complex cultural scenes, the genera- 

'tion^of imagery from a scene constructed according to the DMAAC specifica- 
tioh.s, an assessment of^the resultant diurnal and seasorifil imagery , and the 

/relationship between the'model data requirements and the resources of the. 
data base. 

» * *• 



■ ; • • ■ •• ■ / ".^ ■ : ■ ■ ... ■ : :• . ■■•I ■> 

"^^^ foil owing -sections -01^ the report ."deserve, .-th'-e scene, da'ta'base; 
the tonal, model ^ the simu^l«ted imagery, and. th£. impact on t|ie^ DMAAC data 
. ba^p. The techni,cal. discussion is. follo^ iy t^e conclusions and* research., 
'•'•recommendations. '< • ' \ • .. . " ^ 
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II, DATA BASE GENERATION . . . . . " 

In order to gener^is' LLLTV and .fLI^R imagery, , the construction of a 
data base for a ^eCvted target area according to the DMAAC specif icartioh 
was required^. The data necessary for the simulation of siensor iir/agery 
are divided into three basic typhs. The first type of <^ta consists -of 
a .three-dimensional geometrical description of the features contained in 
a scene. The second type contains descriptions of the physical make-up 
of typical surfaces in the geometric data base and is needed to simulate 
different types of sensors; for example, for the FLIR sensor, the .thermal 
conductivity and capacitance of the material must be known. The third 
.type of data describes environmental conditions for a particular scene * 
^nd: includes temperature of the sun; sky, air, wind speed and di^rection, 
and precipitation. ' 

2,i> Varget area selection 

The first task of this -effort*" invol ved evaluation of the 16 target 
•areas in the AFHRL Seasonal S'^'nsor Handbook* The data for this handbook 
were taken from the tower of the Avionics Laboratory bu-ilding at Wright; 
Patterson Air Force Base and spanned the area between the flight main- 
tenance hangars of Area C to the Trebein radar dish. Scenes in these 
recordings range from nearly all: natural content to a mixture of cultural 
and natural 'features with rather dense cultural content.. The scene 
showihg the hospital was chosen to be the subject of this study because 
of its variety of /cultural content. This scene is shown in Figure r as 
seen from the tower. 

The selection of the hospital scene required -data base model ing of 
a . strip approximately 2-1/2 miles long and 3/4 mile wide, parallel to 
the railroad tracks between Dayton and Xe;nia Avenues in Fairborn, Ohio, • 
as shown in Figure 2. . / 
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Figure. !. Th^e Hospital Scene 



i 



Photographs taken from the tower at the Avionics Laboratory were the 
major source of data used in selecting important features for the construe-" 
tioh -of the. data base. Since not , all details of the stene were vfsible 
from that view, some, of the features were omitted. - Effort, was concentrated^ 
on the haspital (1), officers .family housing" (2), the base heating plant (3), 
the Southwest Portland .Cement, Plant (4), and the Maple Avenue bridge (5). ; 
Exact relative positions, of these features were obtained from a geographical 
• map of Fairborn and a map of the Wright-Patterson. Air. Force Base, Dimensions 
of the features were extrax:te,d from the available blueprints. In addition, 
a number of photographs of the buildings and the surrounding areas wer«e 
taken as an interpretativje aid in .the data base construction. 
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Figure 2. Seogrdphical Map of the Section Selected for the Model 
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A data base was prepared for t'he chosen scene a*t each of three 
levels of detail v These corresponcf to levels lA, IIA.'and 1 1 IB. "of the ' 
1974'DMAAC specification. Level lA has a resolution of 500 feet, level tlA 
has a resolution of- 200 feet, and level IllB has a resolution of 50 feetJ- 
The data bases for tfrese three resolutions were prepared in accordance 
-with guidelines in,"Specifications for Digital Radar Landmass Simulator."; * 

2.2 DATA BASE CONSTRUCTION ' 

^ 2.2.1 Level I A " " -V..- 

•'.j For level lA, only 12 features were modeled and their "feature 
anaTyjsi-s data table" is incTuded in Figure 3. Dimensions of tl^.fe'J / i 
hospital , the steam plant, and' the cement plant fell below the mini nkjm 
requirement of 500 ft heigKt and width for level lA and, therefdre, 
were not modeled to scale. Instead, they were modeled as point 
features, with diameters the same as their convex covers, as suggested'^ 
by the "isolated structures criteria" (100 ft adross and 500 ft d'istant 
from ottter groupings) , • 

Each of the two groups of visible family 'duplexes on the base were 
modeled as homogeneous blocks of uniform height, since the dimensions of 
individual houses fell below the minimum si2e. criterion for level lA, 

All of the features were modeled as three-dimensional polyhedra with 
a base outline as indicated in Figure d. The remaining two features, are 
Route 4 and a ground plane, both of which are modeled as planes. The 
equivalent sensor data base for level lA consists af 111 surfaces con- 
sisting of 188 vertices and 544 edges. 

2,2.2 Level IIA • 

At Teyel IIA, ^21 features met the resolution requirements. Their 
feature analysis data table is shown in Figure 5. Six isolated fea- 
tures were modeled with their basic outlines (as opposed to poijit 
features). These were the officers housing unit modeled to scale. 
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Figure 3. Level lA Feature Analysis Code (FAC) Sheet 
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the lowe revel of the hospital with predominant he-ight of 25 ft, the upper 
level or^^ 41 tne lowert^ level with predominant height of 64 ft, the 
heatingwMik ar^^ the two largest sections of the cement plant (the kiln, \ 
finish T^^^P and rock storage as one and the silos and pack house as 
' "anotherjlllHpe sixth isplated structure is the Maple Avenue overpass* Its ' " 
dimensions (420 by 60 ft) fell within the minimum requirements for level IIA, 
anjl the overpass wat thus modeled to scale. * , - - " 

However, the shape .of the bridge/did not confirm to the DMAAC 

specifications. It' is supported on. both efids by an earth embankment which 

<* ■ " * . .. . . ■ . 

fell below minimum size requirements for level ^IIA^and its shape cannot • , 
be modeled as a set of surfaces paralleT.(or perpendicular) to the ground ' 
plane. This accounts for the "floating ^in' the air" effect the modeled ^ 
bridge has. 

Some difficulties were also encountered when modeling the lower 
portions of the hospital. The height differential of the terrain has ^ 
the effect of an additional floor being exposed on the west side of the 
structure; whereas, it is covered on the east side. t 

Remaining features are homogeneous blocks of residential areas,- 
tree covered areas, a park and an athletic field, two bare ground areas, 
and the base plane. Outlines of all the features in leveKIIA are shown . 
in Figure 6, and the areas between the homogeneoLr blocks give outlines 
of tbe major roads. Altogether, the model at level IIA consisted of 131 
surfaces "Constructed from 611 edges and 285 vert^res. 

2.2.3 Level IIIB / 

Level IIIB contains a data base with the finest level of detail. 
Again, features were- selected from the single view from the tower and 
do not incliide all the 'detail allowed by DMAAC specifications. Despite 
that, there were 87 individuaJ-^atures (some of which had several parts). 
The feature aPialysis sheet is shown in Figure 7. The first 37 features . 
were the family duplexes on the Air Force Base grounds. Most of the 
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Figure 7. Level I I IB FAC Sheet 
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. individual tiouses'Were nat far enough apart to satisfy tfcte DMAAC "' isolated 

structure. (criteria'-', (^o'ft apart) . ' Instead, 'they were .comlained 'intp groups 

running parallal to the. streets^. Th'e' rooftops -of the duplex€^" arp/the 

/ m^jor visible part of .the group- o*f hoiisTng on .the A.ir Fo-rce Base'; ■ \ ^ 

Vjiowever, since the specifications do" not allow for ;s loped roofs, the. 

- model does noj: accurately reflect them.- ' ' ' ' 

Feature 23 is- the "officers quarters" building, which wa's'modeled to' 
. scale. ^ The hospi.tal -(feature 24^) was modeled as fdur separate levels -wit-ti^ 
heights of 9^ ft, 64 ft, 32. ft, and 25 ft, as suggested in' DMAAC'"l^eight 
Continuity Specifications." "Terrain elevation Resented difficulties for 
the bottom lever as they did during the construction at level IIA. 

* I • • (. . . - , , » 

"~ . ' , ■ . r ■ . « « • 

Although the metal shades over the ^outhwesj; walls of ,the hospital ; 
were quite noticeable on the photographs of the scenes, they were hot 
•included in ttie model. Thi^ was, done since the specifications do not 
nave clear instructions for modeling significant details dn sides of 
buildings. All the features carh'only be mofleled as DMAAC "areas wi-th' 
' predominant height." ^ • ' ■ * 

The heating plant, ""service station, and church of Kauffman ^nd Powell 

Avenues were modeled to scale (features 40, 41, 4.2). The Mrfple Street over- 
-. pass was modeled as in lev.eT IIA and with the same 'assumptions (feature 43) .* 

The cement plant wasdivided into 11 individual structures, each of which 
_ had outlines to scale:. (44 - silos, 45 - pack house, 46" - warehouse, 
. 41 - water tower, ;48 - smokestack, 49 - precipitator, 50 - kiln building, 

51 -. finish mills, 5? - rock, storage, 53 - laboratories and offices, ; 

54 - shops amjl stores, and 55 - railroad tracks) . . ' 

• Again, much of the significant detail was not modeled since speci- 
fications do "not allow for sloped roofs or detail on the side? of the 
•Sh^uildings. Furthermore, the rock storage building i-s a str*ucture witti • 
a sloped roof and without one wall; thus , its model does not reflect 
its actual shape. • 



« 

• Refiiaining features are iioHiogenebus blocks'pf either^esidgntiarl / 
, areas or. tree-covered , areas, grass patches, and bare dirt fields. As 
tn , level IIA, the base plane is covered' by the homogeneous area-s "so as to. 
form outlines of the major, streets. Figure '8. 'The level, IffB data baie 
consisted of 684 surfaces maxle but of "3203 ^dges and 1155 vertices^ . . 

2.2rA '.Material Data ftase , - 'i ' '\ ' - ' 

~ ~ ' ' •*■-■' ■ ' ' ^ ^ ^ , *' 

As; stated in the introduction, ttie .approach taken in this study was 

to determine wjj^trinformatiorijs heeded V> simulate the imagery in the 
handbook. This requir'ed 'the derivation of a tonal ' model for ttfe simula- 
tion of .^passi^e sensor imagery'. The tonal model will be described ir> 
S6ction III. What is of concern ^in thi.s section is the additional .daja 
_ba5e descrj^ptors required by the tonal modej . ^ / 

For each of the sensors, var.ious surfaces in the geometrical data 
base will attain different shades depending on th6 physical properties 
of the surfa'ces they represent. ,For the L-LLTV simulation, orlly the 

^surface reflectivity close to the'visiblfe band i-s used^/^9ur4ng the shade* 
assignments For the FLIRynotiel, npit only is the sufface condition / 
iirtportant, but the composition of walls belonging to the scirface is 

^also necessary. This additional i^brmation was added to -.the data ' 
base. A list of available surface material types is shown in Figure 9. 
To completely describe the hospital scene, the material type and reflec- 
tivity of each surface in the scene had to be incorporated Tn the data * 
base. This was done through -a four-digit word. The first two dibits v 
are the material type and the last two digits are the percent reflectivity 

,The additional feature description is given in Figure 10 for each feature 
in the datS base for levels. lA,, HA, and JIIB. for example, the walls of 
the officers housing, feature 23 at level IIIB, has a code 2040, meaning ' 

^they are made of brick, which is 40% reflective" at ttie visible bands. • ' 

The need for these additional data base descriptors and their use' 
will now be described j'n the discussion and derivat'ionf of the thermal 
model . * , * » ' ' ■ ' ^ " * 
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. ' . ' ■ ' MATERIAL TYPES . 

FEAJDRE . FEATURE DfSCRIPTI ON ROOF • WALL • 

NUMBER ■ . . MATERIAL flATERIAL 

■ . " " . LEVEL lA • 



1,6^0 


; .Trees' - " 


3314 


3310 


'2,4,7,8,9 


ResidentiaT Suburban 


3030 


3630 


3 ' . . 


' Commercial 


3430 


3434 


5,9 


Industrial 


3230 


3238 


11 


Roads ' , 


8038 


. 8038 








f 












LEVEL IIA 




f 




Trees^ 


• 9614 . 


i 9610 


2,5,8 

t • . • 


Suburban Residential 


3030 


3030 


3 


Residential Inner City 


' 3130^ 


/ * . 3130' 


4 


Government Buildings 


3430 . - 


3430 


6 


' Industrial (Heating Plant) 


3230 


3245' 


T2 * . 


Industrial (Bridge) 


3230 


3238 


12-14 


Industrial (Cement Plant) 


3230 


' 3238 


15;i6,17,18,19,20 


V Terrain 


3320 ; 


3320 


21 ' ^ 


Rpads 


8038 


8038 
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Figure. 10. (Continued) 
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III. -TONAL MODEL DEVELOPMENT - " ' 

This section presents the electromagnetic and thermal models used 

to produce LLLTV and FLIR sensor imagery. It shows the overall model 

. p ■ , 

- consisting of a scene, environment, and sensor as implemented in a 
software system. # 

3A THERMAL MODELING 

A thermal mod^is the central part of the passive infrared (IR) 
simulation [2]. Thermal sources in the exterior environment and interior 
of objects are considered by this model* Thermal paths from sources to 
surfaces Tink the exterior and interior environments to each surface 
, in the scene. A finite (difference madel of the cross section of each 
* surface (typically a flat slab), with the environment paths as boundary 
conditions, completes, the thermal model, 

3.1.1 Modes of Heat Transfer , 

There are three modes of thenrfal transport. They are conduction, • 
convection, ahd radiation. These are only statisticaJly understood' 
phenomena in that they are the result of a large number of statistically . 
measurable elemental processes: For^present pi>rposes, i$ -i 
only to under&tarid where eich path octfcirs and how to cMracterize-^^^^^ 

Standard equations for each thermal transport method /a j^fe^:^^^ 
the developed thential model. For conduction in' homogeheoiis ^taVfdria^^ 
media, the equation of therrmal transport in a three-dimensional Cartesfan 
coordinate system is 



? 2 2 
3 T ^ O + III = £ 31 (1) 

3?" . 3y^ 3z^ ^ ^* . ' / 



/ 
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where ''K" is a constant thermal conductivity, "c" is ^ constant capacity 
per unit Volume,- "T" is the temperature at; (x,y,z, t) , and "t" is time. 
The thennal model deals only with thermal conduction across a flat-slab 
in one dimension; hence, the conduction equation rieduceis to.J: 



£ il * io\ 
K3t ■ . (2) 



This differential equation lacks boundary conditions. These condition 
are convective air and radiative heat loadings from the exterior and 
interior environments, plus an initial temperature cross section in x. 
Th^s leads to the equations of convective and radiative transport. 

For convective trarvsport across a planar surface at x = 0, the 
equation is . . 



^<T^ - T(0:t)) -c'-^ (3) 



/Where h -is the average convection cpnductance, per unit surface area. 

coefficient based on surface condi_tion$ and / 
l^^^^^^Hv: It: is the cumulative- effect of wind speed and direction 
Wp}0^e.to the swrface orientation, surface^ size, etc. The approach 
;>ped Jigfe is to assign a value based on. ideal, conditions such. as^,laminar 

combi ned. with standard f omul ae using Reynolds and 
liu$seU. numbers. Feature shape, surface type, and orientation, and, .; 
environmental conditions determine these parameters. . ; 

For radiative transport across a planar surface at x = 0, the equa- 
tion is 

...... ■ • , -s 

N 



E (a, E,) - CT^ = c i^ti* W 



where^there are '"N" sources of Trradiant piwer "E." r "€j" is the average 
thermal OTtttance, "a.'V is the average absorption for irradiant energy of 
"i's^* spedtniim, and "a" is the Stefan-Boltzmann constant; 

3.1.^ The Viewing Factor . ' ^ • 

^ The characterization of radiant paths requires a measure of energy 
exchange: the viewing 'factor. The viewing factor is used when two sur- 
faces radiate to each other and one wants to calculate the amount of 
energy ^interchanged.' The amount of jBriergy reaching surface 2 f rom ^ 
surface 1 can be calculated as follows. ' r 

Surface iS characterized by it^ diffuse radiance. For each dif- 
ferential element dSj, a portion of the ra(^iant ielri^rgy hits surface S^. 
This part of the radiant energy of surface Sj^ has to be surraned in order 
to compute the radiant interchange between them. This suiranation can be 
characterized by the equation . 



1 / / cos cos 69 

'uz-li r \2 

I . Trr 



dA. dAg (5) 

7rr~ 



where "9:" is the angle between the nonfial to dA^- and the l*fne segment 
connection dA^ and dAg* "r" is^the length of the line segment connecting 
dAj and .dAg, and ^A^-" is the area of surface S^. . The shaipe factor is 
used in radiant interchanges where the surfaces have constant radiance 
across themselves as .viewed by each other. In -such a case, where the 
radiance of surface is and the radiance of surface is In 
(as viewed by, each other), the shape factof is used to calculate* the 
net radiant interchange b/ / 
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For thennal modeling, the shape factor can be used with surface 
temperature to compute heat^ loading power by 



The above procedure demonstrates how radiant paths are modeled in the 
finite differende thermal model. 

3.1.3 Environment Paths 

A simplified environmental model is used. Figure 11 presents a 
schanatic of the model. The ground plane, sky hemisphere, sun, and the 
four diffuse sotar quadrants are the exterior radiapt paths to a surface. 
The- atmpsphere is the exterior convective path to a surface. The inter- ^ 
iqr side of a surface has a free convection path to the interior air and 
a radiation path. ' ^ 

The ground plane is a horizontal plane of infinite extent with con- 
stant diffuse thermal radiance. It is a grey body with a spectral peak 
near 10 microns. Its diffuse radiance is related to temperature by 
Stefan-Boltanann's Law. The ground gan be dry, wet, ^ or snow covered. 
Under wet and snow covered conqitions, it closely approximates a black 
body, ^so the ground condition dbe'^s not have mijor impjact on its thermal 
emission. When the ground is snow covered, though, its maximum temper- 
ature is 32*'F* 

s 

-"—The sky is a hemisphere over the ground plane. It is a close approx- 
imation to a black body except in spectral bands where the isky is highly 
transparent. It has constant radiance over the hemisphere because almost 
all the received thermal energy is from the bottom 1 kilometer of the 
sky. The sky's spectral peak is near 10 microns. The presence of clouds, 
especially low, rain-bearing clouds, increases the sky's radiamt, or 

r . 29 ^ 
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.Figure 11. Schematic Representation of Exterior Environment 
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" apparent; temperature" because cTbu^ those bands where the 

sky is highly transparent, .the sky and ground have ^ilar thermal 
emission bands since their temperatures are similar. 

The sun is a*" SDOO^'K black bocly. It is a source of both direct and 
diffuse radiation. The direct radiance J s modeled as Wing emitted by a 
point source, while the di ffuserradiance ji$ tiwdeled as coming frum fbur 
sky quadrants. Each quadrant. is a constant radiance : diffuse surfax^^ 
approximates the diffuse solar radiation falling from that ^quadrant of the 
sky. The spectral distribution of solar r^idiation is significantly dif- 
ferent from that of ground and, sky. About- 99 rpercent of solar radiation 
is between 0.2 and 3.0 microns which contains the visible wavelengths. 
However, the bulk of the received'^adiation from the ground and sky is 
at wavelengths greater than 3 microns. 

The atmQsphere is modeled as ,a gas* Its temperature and wind vector 
determine the rate of convection between a surface's exterix)r and the 
atmosphere. Reynolds and Prandtl numbers are ^sed to estimate convpction 
rate. Precipitation influences exterior fieal; transfer since Incident Vater 
is heated bj^ conduction to the surface temperature. Based on precipitation 
ratej wind vector, and surface orientation, a certjiin quantity of water 
hits a square meter of surface-Per hour: ^fhfis Water is litea ted from its 
perci pi tati on temperature to the^^rfjce tenip^ratu^ r^suli^ng 
heift transfer and evaporation. " ' ! ^ 

' The interior side of a surface has^a free convection path to the 
interior: air. The Interior air may be air conditioned {i.e., iTi the 
ramge-eS" to 75''F), heated (I.e.', keipt-aii6\/6 eS-Py, urrheated (i .e. , 
enclosetd but not heated or coo1^dy,"^oKop^ri to th'e'-^crufside. • Based on 
the condition of the fntefibr atK^fr^^'t'onve'tlion from the'iriterior 
surface to the enclosed air mass -can occiir. " ^ 



, 3.1^ fhe Finfte Difference Model 

The geometric description of objects for the theririal model consists, 
offset of points, line segments, and planar polygonal -bounded surfaces, 
^h is modeled as either a flat layered slab, a cylinder, or a sphere. 
Periodic layers like beams or studs can be modeled. Layers are sectioned 
with thermal capacity within the section and thermal- resistance to the 
adjacent sections ^issigned. The edges of. the slab form the boundary of 
the slab with its exterior and interior environments. Boundaries or 
edge effects between different wall cross sections are ignored- This 
is the basis of the finite difference model. , . , • 

•- / ' 

The finite difference method is a powerful' method "of ap'pijoximating 
rtany types of partial differential equations. Applied to thermal model- 
ing, the method permits quantitative computer modeling of the walls, 
roof and ce'iling*, and also interidr and exterior convective and radia- 
tive fnteractions of- a target feature. 

Figure 12 shows three typical "wall cross sections: uninsulated 
wood, insulated wood, and brick. The wall models interface with the 
environment at their interior and exterior, where convective and radia- 
tive transfer occur. The interior wall models include conduction through 
s^olids and' convection arid , radiation acros^Nair gaR^, the para IT el ther- 
mal paths of wood studs and either air gaps -or insulation are modeled as 
parallel one-dimensior|al paths without cross paths included. This has 
been found satisfactory for mea"sured cases. .' 

The finite method hats a theoretical foundation in approximation 
theory using function spaces. The theory of application has been 'devel- 
oped and includes criteria for a well -formulated model in terms of con- 
sistency, convergence-, and stability. A finite difference equation 
(F.D-E.) can be written for a node "i" connected to two other nodes 
in the model : . ' . 
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Consistency demands that as the incrraTents of time and distance approach ^ 
zero, the f/D.E. approaches the partial differential equation. Rearrang- 
ing (9), we see that^ 

(T.^j(N)-T^(N)J - (T.(N)-T. (H)) ' T.(N+r)-T,( 

• „ • . . ■ 

As AX,-At 0, this approaches the differential equation for conduction, / 
As contrasted with- consistency, which demands that the difference equation 
at each nodp approach the partial .differential equation, convergence demands 
that the solution of the difference equation at each node approach the solu- 
tion to the partial differential equation at a/. At ■* 0. For most' cases 
of thermal modeling, it is sufficient to prove the third criterion, sta- 
bility, since convergence of solutions, is then assured. ; Stability refers 
to the behavior of the numerical solution calculated ',by finite wopd 
length computer. The stahili-ty criterion gtiarant;ees that the truncation-^ 
erroVs of \iterating the F.D.E. do .not propagate. - , / ' ' - 

Using the model, the temperature at discrete intervals within 
a wall. and the heat flow from (or into) the exterior surfaces were 
computed and are illustrated in Figures 13 and 14.^ v 

Figure 13 presents the temperature of each node of a wall a§ a func- 
tion of time. Tor contrasting thermal behaviors^ an uninsulated wood wallt^ 
an insulated wood wall, and a brick wal^l are shown. Surface materials arid 
their cross sections have a major impact on exterior surface temperature. 

Figure 14 presents the thermal loading of three different environ- 
merits as a function of time for some .standard surface orientations. The 
environment, time, and surface orientation all have significant impact 
on thermal exchange at the surface. ^ 
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Thermal modeling proceeds by first modeling the surfaces* of a target 
feature and then connecting -these elements to the thermal , paths available. 
In thts* approach, walls, ceiling aftd roof, floor, and internal air condi- 



tions arjB. modeled individually. For each feature, the appr6priate standard 
•* * • » ' 

elements are cbnnected by- thermal paths forming an overall thermal network. 

The full method of constructing finite^ifference:^network models of 
target features requires a scheme involving standard component parts with* ' 
indiviflual submodels^. These components arje ,define<^ later in architectural 
classificatiens. ' — ^ j 

3.^ EMITTANCE. ANI^ ABSORBANCE MODELING I 

The amount and spectrum of radiant^ energy emitted at each angl^ by ^ 
a surface is ctoracterized by the surface temperature and spectral, 
directional emittance. The blackbody spectral emissive power function 
af Plank is used to calculate the maximum emissive power as a function 
of wavelength and then •multiplied by the spectral, directional emittance 
at each^ wavelength. \ Plank's equation ilr ^ 



X5(e2 



1 

c;7xr 



(11) 



-1) 



where "E." is monochromatic emissive power in (W/hr m y). 



X" is wave- 
is 0.69544 X 



length in microns, "T" is temperature in Tcelvins, '"c, 
12 2 

10" ■ (W cm )i and "c^" is 1.4388 (K cm). The approach used here is 
to estimate an average of the spectral, directional emittance oVer a < 
-band, weighted by the expected amount of energy output at each frequency 
across the band. Then, a diffuse emissivity is estimated from the normally 
directed emissivity. Emittance is th^ used in either thermal modeling of 
radiation paths or prediction of tfiermal infrared signatures. • 
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Emittance is a statistical property not tompletely determined by * 
the surface matew^%^'^^^^^ including roughness and degree' 

of oxidation, have a major impact on emittance. * The volume of material 
just under th.e' surface also contributes to emittance in pohcbnductors. 
Figure 15 shows [some, dir^ctfanal emittance curves of prepared clean 
surfaces. The cjurves approximate diffuse emissive sources except near 
grazing.ahgles, ,: In , addition to the amount of energy emitted, there can 
be*^ degree of ijolarization to the energy, especially in the millimeter 
bandr • ' • •■ ; ~ ■ . • ' / 

♦ 0 /The basis of thermal emission modeling is Plank's law of spectral 
blackbody emission. From this equation several radiation laws can be ■' ■ 
derived. Integrating Plank's equation, one, finds .the total blackbody 
thermal hemispherical emission is 

• \ • . ' . # ■ ■ ■ 

where "a'; is the^Stefan-BoltSmann constant. Wien's displacement law can 
be derived -from Plank's equation. It is / 



where "C3" -is 0.28978 (K.cm). < - . • - 

At the temperatures of most objects on the ground, the bulk of 
emitted energy is longer than 6p. There is a need to model energy* 
reflected and absorbed ^n the emission band of ground features (i.e., 
roughly 3 to 15p) in order to compute radiant pfiths for thermaV modeling. 

In addfti^n to the 0.2 to 3. Op and 3.0 to 15. Op bands, there is a need 
to model emittance, absorbahce, and reflectance for 0.8 to 1 .Ip' (LLLTV) and 
8 to 14p (FLIR) sensor band§. Emittance, absorbance, and reflectance func- 
tions are required for eatfh spectral band of interest, i.e., those bands 
which are significant indescribing the ground, sun, or a specific sensor. 

• ■ ■ ■ •' 
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Finally, Kirchoff's Law plays an tiflfortant role here ip relating 
• emittance (e), absorbance (a), arid refTectance (p) functions. k?rch6ff' s 
Law equates emittance wHh absprb^^ oi^ fqui1 il)r^tum. 

Thus, under. the condition that fransmittance is zero, a =^ e = l"^- and 
aW three radiation interactions are related in eacft band.,' r 

3.3 "REFUECTANCErWDELING 

The amount of reflected radiant enjergy in any direction as a function 
of the incidence angle of the irradiant energy is ' characterized ^^^^^^^^ 
direGtlQnal reflectsrice: : The btdirectiorial reflectance is the; ratio of 
thevmcremeata^T radiance dlyXe. ^&.) over the incremerital irradianCe- 
dE^CO^'.- ■ .v'-:!-c i.: r;. ...v, ? .:■ 
■ '■■■^-'c-v . :y - ..: J'- A- ' ' ■' ' -> ' •■ ; ■. 

The bidirectional reflectance is a statistical phenomenon. Most real 
surfaces 'are too complicated to have; a single constant reflectance. In 
these cases, the bidirectional reflectance represents only the mean reflec- 
-ta'ntS€f for the givert iYicidenc€ and ;refTectance angles. Since there *s a 
signifldsnt refTectance variation ■ point-to^potntv^along tfje^ surface, 
care ts pjBCfessairy when Interpreting mG^asure^ bidirectidnal reflectances. 
One a>prtfa?(^ is tdi^fl Iter meiasured data to eTimlnate high frequency 
angular* changes -i n ^'^jef Tectarrce si nee theseVcontai n much of the vari abl e- 
oyer-^the^isurFacfe li^T/onnatlon (i^e*. ,'^rioisfe)>; ^ : ■ 

Polanzation is an important factor in reflectance^behavfbrr Polar- 
izations parallel and perpendicular to the plane of incidence have signif- 
icantly di f f erent ref 1 eci:f on behaviors so tftat aif ter fef 1 ectlonsv 1 iglrt 
originally unpoTarized can gain a degree of polarization. 




The simplest model of reflection is Lambertian (i.e., cosine law) 
diffuse reflection. This is a constant radiance reflection over the 
hemisphere above the surface;, Many simulations use this simple diffuse 
model since it appears to be an adequate description of the scattering 
mechanism in many cases. 

■ The next step, which was not within the scope of this study, is to 
develop a capability to include specular reflection in the model.. 

- 3.4 ATMOSPHERIC MODELING 

The atmosphere interactswith radiation via scattering, absorption, 
and defocusing. In the model developed for this effort, defocusing has 
been ignored, while scattering and c^sorption interactions are easy to 
characterize by the attenuation coefficient a. When applied to a radiant 
surface' at a range, r,, the transmittarj^e in percent is e"""^. If the sur- 
face has radiance:! and the atmosphere has apparent radiance I at in- 
finite range,, the observed radiance is 

^0 = ^s^""' ^a (1 - • (15) 

This is simplifying the underlying phenomena. In-order to be accurate, 
attenuation should be a))pli^d to radiances wavelength by wavelength to 
account for spectral inhomogeneities in the atmosphere. Yet, that would 
require a prohibitive amount of computation. When atmospheric attenua- 
tion, is involved, radiant energy is computed for spectral windows, and 
an average attenuation rate computed over^eacji window, wihen the paths ' - 
are sfton, as ^betweeVi nearby surfaces, .atmospheric attenuation is ignored. 

3.5 A METHOD OF CHARACTERIZING ARCHITECTURAL. SYSTEMS 

In order to perform thermal modeling, cross sections of walls and 
roofs and , interior air volume are needed. These components provide 
thermal -memory., t^Wertunately, it is impractical -to determine the 
cross section of each surface in a scene. 



The approach taken here is to construct standard cross sections that - 
are labeled and can be entered into a.model as required'. The valid-ity of 
thfs approach is based on all the j2onstra>nts placed on building designers 
(and nature). Building codes, architectural methods taught in school, 
available construction materials, and environmental conditions all limit 
the types of design. 

The first step is to characterise a feature's overall construction 
type. Figure 16 illustrates the top level of a construction classifica- 
tion tree for buildings. Other classification trees can be constructed 
for pavements, natural areas, etc. For buildings, the overall design is 
typically easy to obs^ej*ve from photography or i? known via scene mensura- 
tion. Other types of itonstructi on, such.as pavement- or industrial 
processing structures, are also reasonably easy to visually identify 
and subsequently model thermally. ' 

' There are many subtrees of this classification scheme. Any building 
can be classified as either a single category or a mixture of categories, ' 
then each 'identified building component (iVe., walls, roofs, etc.) of the 
geometrical data base can be thermally modeled by a standard cross section. 
In this way, the total building is thermally modeled in terms of its .con- 
struction, interior environment, and componenV surfaces. 

It is possible to construct a thermal model of a feature by (1). clas- 
sifying it in terms of standard categories, and (2) applying ideal material 
cross sections (corresponding- to the named classification) to identified 
feature components (i.e., wall, windows, and roofs) in the geometrical 
data. 

3.6 IMPLEMENTATION OF THE TONAL MODEL 

The thermal model is performed in two processes. In these pi^ocesses, 
the geometrical and thermal properties of individual objects in the scene 
are processed together with the dynamically changing environmental con- 
ditions in order to assign, shade values to the external surfaces of* the 
oyects; The entire flow, diagram o^the thermal model implemisntation . 
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isshown in Figure 17. I it the first process", a'matrix of coefficients 
is tonstructed with an entry for each wafv type in the .scene. Then in 
the second process, each of the surfaces An^.he. scene is assigned a shade, 
using the above constructed matrix. > ' k ' 

3.6.1 Construction of a Matrix of Coefficients 

For each of the availabTe wall types (Figure 9) aifid, each of the 
standard orient;atfons, a set of four coefficients is^ generated,. To ob- 
tain each set of the coefficients, a temperature vector. isj computed for 
five typical ' thermal reflectances of^he; current wa IT type. ; The first 
element of each-of the five vectors is used to find a' foucth tiegree 
polynomial by a least squares fit raqjthod. Coefficients of. this poly- 
nomial are the next entry in the matrix of coeff iciehts." " - 

To allow the tOTperature of *lie ^all ^ 
is necessary to ddnsidisr riot onl^j^he current cc^dftidri of the environ- 
ment, but alsojthe conditions of the surrounding environmeafj during 
several preceding hours. Thus, the tefniperature vector'.is computed in^ 
several iterations over the time increment.. 

The temperature vector is set initially to have -its first element 
equal to temperature of the interior air and is obtained-from the "wall 
type file." The last element of the.. temperature yec^ initially 
set equal to the current temperature of exi^erior air wh-fct] is obtained 
from the "environment condition filie,'- and the rest of the. elements are 
linearly distributed between the fi"rst and last elements. Each of the ^ 
elements corresponds to. a "node" of the current wall type/ The heat 
flow through the wall depends on the thermal properties of. the wall 
and the external loadiiig and can be ana^tically represented by a 
partial differential' equation. Furthermore^ using the first few - ^• 
terms of a Taylor polynomial for two variables, the differential 
equation for temperature (time, position) can be reduced to a set 
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of linear equations involving, the thermal capacitance, ronductance, and 
'.time increment as coefficients and temperatures at various nodes as 

variables. This determines a matrix, and the system of linear eijuations 
■ is-. solyed by applying its inverse to the initial temperature vector, id 

inclijde^the influence of the change of the environment on the exterior 

node (surface), an extra load node is added. The load depends on:. 



(1) Sun position with respect to exterior surface of the wall 

(2) ' Amount of s*:y "visible*" to th^ surface. i^ 

(3) Amount of the ground "visible" to the surface. 

(4) Air temperature, .wind velocity, and precipitation. 



To summarize, at each time ;^iii±eryal , a temperature vector is fouiid where 
loading due to the sun, sky^'and ground radiance^ and tKe air convection 
is added as an extra node. Such a vector is generated for each of the 
typical thermal reflectances of the current ^all type, and the last ele- 
ment of these vectors is used to determine the coefficients of the • 
fourth degree polynomial. This is^done for several orientations for- 
each of' the wall types. Th6 wall type and orientation are then used to 
index the matrix of coefficients. 

3.6.2 Shade Assignment 4 

Once the matrix of coefficients is constructed for a specified time 
of the day, the actual assignment of the Shades^to the individual surfaces 
can be performed. The wall type, the normal vector, and the thermal reflec- 
tance ar^ serially retrieved from tfie geometrical data base. For each|of 
,these, the three closest standard orieatations are found, since the matrix of 
coefficients contained only a fixed number drjentations. Next, coeffi- 
cients for these orientations are used to interpolate the approximate 
coefficients for the gfven orientation. Finally, a fourth degree poly- 
nomial with the approximated coefficients is evaluated at a given reflec- 
tance. This value is then properly scaled and used as the shade for the 
surface. ^ * ' 
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Xhqc:scene, daU bases fori the three different 1 evel s of detail Wd 
j*iei|.. construction were diescrtbed^ilp Section 1 1 ,s and the vi sua! aild thermal 
radeiy-were developedt. in Seetion llJ.. In-.^his ,t4ctionV i.magery :generated 
from the data -b?ise and the thermal /visible model s^and the im'agir^y froin 
the actual sensors are presented. • Several sensitivitiesi were selected 
to lllustra'^'e the effect of the changing cdfiditi on of the environment 
on the imagery. For each of the %d Senso/s, LLLTrlnd fLlW, three ■ 
signi f icantly di f fereni''days were cl|bth to demon stra>b" the 'Infl uence .>.,' 
of changing seasqni.' The environmengn conditions corresponding to 
the 3 days are listed in Table 1. The first of the 3 daysp's a cool . 
clear day during the winter, with the direct sun^ radiances, the 
major. component of the. iHumin^tioti and the thermal loading. . The second 
day is a hot cloudy day dMring.,the summer, and the diffuse sun radiance 
(iUumination by the sky) and tliWlr temperature become more signifij 
cant than on the clear, day. The..third day, a fall day, is ah average 
warm day with partial cloud cover. The winter, day produced the best con- 
trast in the images. Four different times during this day were selected' 
to illustrate the sensitivity of the scene to the diurnal variation. 
For this purpose images at 10 AM,"" 1 PW, and 4 PM are included 

for'both sensors ateach of the thr^e levels bf the detail. 

1 Ftgure 18 shows thp imagery taken from the AFHRL -Sensor Handboql^ 
which ^are photographs of the.sq.end taken by the actual sensors. Jn : 
$e.ction 4.2, imagery simulating the.FLIR sensor is presented; imagery 
simuUttng LLLTV sensorS'_is presented in Section .4.3. These two sections 
demonstrate tfie seasonal J diurnal, and level of , detail -sensitivities. 



TABLE 1. ENVIRONMENTAL CONDITIONS ^ 



Season 




Day 


Cloud 


Temperature 


Wind* 


Winter 




7 February 1976 


Clear 


Cold 


Wi ndy 


Summer - 




3 September 1976 
7 October 1976 
— If 


Clear 
Cloudy 


Warm 
Cool J . 


Still 
Still - 



4.1 ACTUAL VS. SIMULATED IMAGERY 



Figure 18 shows the imagery of the hospital scene taken by the actual 
LLLTV (18a, 18c, and 18e) and the FLIR (18b, 18d, and 18f) sensors. In 
the first photo taken by the LLLTV sensor on a clear cool day , ,.the position 
of the sun is clearly indicated by the shades of the hospital wallsfacing 
towards the observer (FiguK;^18a). Notice that the sirpulated imagery 
for winter provides a good approximation of ttje image as seen by the real 
sensor. On the same day,- in the FLIR^imagery (Figure 18b), all walls are 
■qfuite bright but of uniform shade. . Again, the simulated image has the 
same characteristics (Figure 19c). Images in. Figures 18c and .18d are 
shown for 2 PM on a^hot overcast day. These conditions correspond to 
those of the simulated summer images (Figure 19d). The last^wo images" 
are at noon on ,a partially clear warm day which corresponds to the 
conditions of the fall images (Figure 19c). Notice that the resolution . 
of the real sensors far exceeds that of even level IIIB (which is .about 
50 ft). > ■ . > » 

4.2 SIMULATED FLIR IMAGERY — I ^ •• 

4.2.1 Seasonal Variatio n ' - ' 

Figure 19 shows three jmages of level riIB~%t various times of the 
year. Environment conditions dei^ined for t?hj simulation of these images 
were the same as in Fairborn, Ohio at 10 AM, jon 7 February, 2 September, 
and 7 October 1976, respectively. J I " : 



J 



W1 nter - the w^nieri^ day was -a c1 ear cool day witfj' sky and air tempera- 

tiir^s Clojs#?ta 40"^^^^^^ ifadiance wil' tte most significant component 

■■■■■ ''I'Mi^^e:,^^^^^-^ ■ " - ' 



of tH^eat-toadfiiT^^ resultecl In higher sensitivity of theVwalls 

to the <posHipn of the. s.un thil As a consequence, 

valls-^facirig' fn tfi^ directionrof the sun are significantly warmer (iridf- 
catedj&y^a ligli^^ter shade) than are the walls facing away from the sun. For 
example, the sou^h-fac+iig wait of the hospital is quite Warm, whereas the' 
west ^facing wall of \the hospital , which is not irradiated directly by the 
sun, is :iti 11. quite pold (Figure 09a). Also notice that the southwest- 
facing walls of the c^^leni^ plapt, -eyen thougpfethey are not directly irra-^ 
diated, display a. noticM temj^brature gradient. This illustrates the 
conductivity of the Vari^ous materiafs. The. left-most Wall is concrete- ' , 
block; and thus", acts^ a? a fairly good insulator between the cold exter-^ 
lor and the warm (heated) interior. This accounts for its dark shade. 
On the other hand, the right-most wall i^ a metal structure and is a 
good conductor of heat. Thus, heat is being lost t)tmjgh this wall into 
the cold environment, and this is indicated by the Ijg^t shade of this 
met^l wall. 

Simmer - The' en v^ironmental- conditions for the typical summer day 
areJIot and cloudy. The temperature at 10 AM was around 65'', Due to the 
clouV cover, the^alls were not sensitive to the position of the sun as 
tjjey 'were on a 61ear day in winter. All walls- have fairly uniform 
shades since tlie'ait* and the sky. temperature's are the dominating heat 
load components. Sky radiance due to diffuse sun load is higher -than - 
J* n winter, and the roofs of the buildings have significantly highe/ 
temperature than walls. Thus, the temperature difference between thei walls 
(and rt^oj" is tiigher- in summer than in winter. ^Oreover, the overall Jem- • 
peratLTes are? higher than in either falT or winter, as indicated by the 
b/ighter image (Figure igb). 
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Fan The day^ for fall was a reratively clear, mrro di^y. The 
temperature of the air was slightly lower. than in surrtner. ; This causes., 
the walls to^be generally cooler (darker shade) than in the summer 
image. However, temperatures are still warmer than on the winter day 
(Figure 19c). . 

4.2>2 'Diurnal Variation - 

Figure ^ gives the outlines of the major features/ in the scene together 
with "positions of the sun at the^various times of ylay. ' Figure 21 contains 
the simulated imagery of the scene.at 7 AM, 10 AM, 1 PM, and-5 PM, , 
respectively. Environment conditions at these times are'' the same as on 
7 February .1976 in Fairborn, Ohio./ 

Early morning - At 7 Af>l,;the surf is still below the horizon; thus :,r^ 
the only thermal loading iis due to the sky»radiance and. air temperature. ^[ : 
All the walls have approximately the same shades since^ they are all uni- ' 
formly loaded by the air temperature (Figure 2Ta). Notice that the roofs 
of the foreground buildings are significantly cooler; (darker}^ than the 
corresponding walls. This is caused by the low temperature of the sky 
whick forces this roofs to radiate heat'. ^ - " . , 

Mid-morning - By 10 AM, tfhff^un has risen 'above the horizon and is 
the major component of the heat load. This is reflected in the high 
contrast of the Figure 2^b.^ Shades, of the waTrls, especially those of ^ ^ 
the hospital, clearly indicate the position of the sun. South-facing 
walls of the hospital have signifiGantly higher temper'^tiires than the 
west-facing walls. This is because the sun is illuminating south- 
facing wal K but not^the west-faciha wal Is (Figure 20). The roofs ^ 
of the building in the foreg^round a^e at this time also illuminated 
by the sun, causing ci reversal of relative temperatures of walls arfd 
roofs. ./ 

Afternoon - At 1 PM, the sun rises close to its highest position 
(apex). and subtends largec angles with the normals to the walls of the 
hospjtal (Figur^e 20). This results in'brighter ground surfaces and 
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roofs,- and insignificant change in the shades of the hospital walls as 
shown in Figure 21c. However, where at 10 AM the walls of the cement 
plant were not directly illuminated by the sun, they. do face the sun, by 
.1 PM; This accounts for the lighter shade of the cement plant, wal Is when 
compared tevthe shades at 10 AM. Also, notice.that the metaJ^wa-11 (right- 
n0t) is wartmer then the concrete block wall .(left-most). This is due 
to h'igheV* thermal ■\^nducti'vity of the met'aT when compared to the conduc- 



tivity|pf the concrete, ^ • 

Evening ^- 5 PM, the thermal loadings due to the air temperature 
and the dffect sun radiancy are about equal This produces the "washout" 
effect in Figure 21d, . since the waTls are not tfs sensitive to the posi- . 
t,iQn of>the sun as they were in ,the middle of the day. Generally,' all 
the sha^esare .significantly 1 ighter than during the morning; hours , 
since they had time to warm up throiip^the, day and the thermal capaci- 
tance aliows'lthem to, retai^'the heat into the early evening houV 
Also, notice the shade reversal , especialljr on the heating plant walls. 
■When. compared to the shades thr-ouglt^the rest of the day, walls facing . . 
south and westv reverse their.shades dur.ing early evening. This is 
caused by the passage of .the sun throughout/the 'day as indicated in. 
Figure 20. ' • - - v. ' ' \ 

4.2.3 Levels lA arid IIA r .. >. . ^ - \ 

Images at' the level lA have nb lairge siir/aces facing the obseryer; 
thus, during the^'diurnal vaniatioriV the sufi position cannot be detected 
as easijy as at the ^^her levels of detai|l. Only the overall change 
of th$ er^ti^onment throughout: the day can be trac'ed.in these "images. ■ 
Sarfacels '.jn the scene get, progressively warmer as ttfe'sun rises and 

warms u'p%he suitrounding air as indicated by the gradual brightening 

>' • . ^ ■ ... ■ ■ . ■ , ^ ■ ■ ■ ■ ' _ • 

of the images (figure 22). . . ' " ' 
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'■At level IIA, there are more individual features and both the -^position 
of the sun and the temperature- of the environment are quite apparent* ^ 
. (Figure 23). SQutb-facing walls of alT three jnajor individual features 
(hospital," heating plant, and the cement'plant) are warmer tha'h the walls 
facing west, and the images become gradually lighter as the day. progresses , 
just as they do in levels lA and'lIIB. 

4.3 SlflULATEP'LLLTV IMAGERY - ' - 

Only two parameters are used in th§ visibly model for the computation 

of the surface shades." The more significant of the two is the direct sun ' 

radiance." The other parameter is' the diffuse sk>i radi-ance which is a 

magnitude smaller than the direct, radiance. SUrfices in the mod^l were 

thus very sensitive to the sun positfon, and all the images hav€ in -^- '^"^ 

general higher contrast than the images in the ^thermal simulation. ^ " 

. ■ ■ . " J 

4.3.1 Seasonal Variation ^ ' * 

^Figure 24 contains images at three seasons of the year. The first 
one is on 7 February which is a clear cold dfiy with significant sun rad- 
iance; the second image is on a cloudy day on 3 September, and the third 
daj^is a relatively clear day on ? October. ■ . 

Winter - The most sigr^ificant component of the illumination is the 
sun radiance; thus, the image has high contrast. Surfaces facing away^ 
from the direction of the^ sun are significantly darker than the surfaces 
facing toward the-#sun (Figure 24a). For example, the li^pital walls 
spfacihg south (Figure 20 )^^|r^ighter than s,urf aces facing west. 

Summer - On a cloudy day i Figure 24b, the scene is not as sensitive ■ 
to the position of thp suis^s on a-clear day. This is because the 
direct sun radiance is no longer significantly larger than the diffuse . 
sky>raKliance. Illumination thus comes from both of these components, 
and the image^contrast is reduceil Notice- also that the overall 



illumination of the scene is reduced, and the image appears to be darker 
than on the clear day in winter. 

Fall - Figure 24c is the image for a parti ally cloudy day, and thus 
the direqt sun radiance Is higher than in the summer cldudy image fut 
lower than in the clear winter image. Again, as in the siimmer scene, 
clouds -in the sky diffuse some of the direct- sunlight and cause the 
sky to appear brighter than on the clear wintj^r day. This accounts 
for the low contrast of the image. V " . / 

4>,3.2 Diurnal Variation 

Images for 7 AM, 10 AM; 1 PM,.and 4 IPM, respectively, are sha^n in 
Figure 25 and\^re based on ,the level I I IB data tase model. The sun 
position and intensity $it all four times correspond to those on - 
7 February* 1976 in Farirborn, Ohio. 



the oaly illumjoation is due to the sky radiance, i 
very similar -shades causing the "washout" effect in 



Morning - At 7 AM, the sun has not yet risen above the horizon,, and 

All the surfaces have 
Figure 25a. At 10 AM, 
the sun is still relatively low and to the east. The position of the sua 
is clearly indicated by the shades of the wall of the Upper sectioji of the 
hospital. Surfaces' facing east are substantially brighter than the' 
surfaces, facing we§t (Rigures 20 and 25b). .However;, the southwest walls . 
of the cement plant f-ace too far toward the weSi and thus are not illu- 
minated. * 

^ Midday At 1 PM, the sun is at its highest point and illuminates the 
scene w'ith the maximiAn intensity^ (Figure 25c). South-facing walls of the 
hospital reach their lightest color of the day as^ do all. surfaces parallel 
to the ground plane. At this time the sun has also advanced far enough to 
illuminate the southwest- facing walls of the. cemept plant and the west- 
facing walls of the hospital (Figure 20). ^ * - 



Evening - At 4 PM, the sun has moved to the west, and the south- 
west-facing walls of the cement plant are almost perpendicular* to the 
direction of the sun'. Thus, their sh,ade reaches the brightest point 
of the day.. The sun has advanced far enough ^ the west that fhe south- 
facing walls of the hospital are illuminated 'ss than the walls facing 
west (Figure 25d). This accounts for the shade reversal of the hospital 
/iSalls at 10 AM and 4 PM/ " 

3 Levels I A and HA ' 




was the case with the IR sensor, images of the leveV lA do not 
^provide much 'information about the intensity or the position «f the sun 
(Figure 26). Only the large surfaces parallel to. the ground plane proV 
vide some' indication of w^tere the sun is.^V. in MrticuU sun rises 

and the angle of incidence increases, surfaces p.arallel t.o the ground 
become lighter. V > • • 

At level IIA,''Figure 27, there are more individual features, similar 
to the IR images'. Howeiver, the position of the^sun cannot be determined/ . 
from the individual surfaces. ' / . - ' , - ^ ° 
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\ . V. IMPACT ON THE DMAAt DATA BASE 

•At this ppint, the major aspects of the project including the data 
base, tonal models, and imagery riave been presented. The intent of this 
work was to evaluate data base descriptors and algorithms for LLLtV^and 
FLIR simulation. In this -section, the descriptors needed in LLLTV and 
FLIR simulations are presented, the descriptors available in the DMAAC" data 
base are discussed/ and a resolution of any discrepancy is attempted. 

5.1 REQUIRED PARAMETERS >' ^ 

: The data required; to run a sensor simulation can be divided into 
three categories: scene data, environment data, and sensor data. )The 
scene data consist of the shape and material description of the sgene. 
The environment data^ consist of the electro-magnetic and thermal situa- 
tion surrounding the scene a[t each instant. The sensor data consist of- 
the performance behavior of a sensor in imaging the scene as illuminated: 
by the^nviroriment. these data actually exist as parameters, i^.e., posi- 
tion, material code,*:radiance, field j)f view, etc. : 

'Parameters used in the simulation depend on the model and, algorithms 
employed and are highly interrelated. As such, it isjmpossible to ex- 
plicitly state a set of parameters that must be in a scene, environment, 
or sensor data base. Some parameters are derivable from other parameters, 
so that the choice Qf parameters used is sometimes arbitrary. With this 
in mind, m;(njmal parameter sets and information derivable from data base 
parameters are useful ways of describing requirement's on the data base. 

Parameters are correlated with many aspects of a feature,' and the , 
utilization of these correlations is necessary in order . to use: the DMAAC 
.data base* Surface and feature parameters are correlated with feature 
use, materials, location'l^ size, age, and architecture. Architecture ' f 
seems to be the best basis for inferring sensor-related parameters. 



It is highly correlated to physical condition (i .e. , accurate, highly 
descript^e, detailed, and compact). I* is simply and accurately 
extracted .from most available intelligence, • < 

5.1.1, FLIR ' " /IV 

The most cjifficult simulation is passive thermal. The sc^ne ligna- 
ture -has a thecal memo'^y of up to several hours previous to'the imag^; 
'thus, the imi^is a result of numerous P^-evdous e\^ts. Figure 28 pre^ 
sents the significant categories of data and expected sources^ fJir ' ' 
simulation places heavy demands on the DMAAC data base and requires a 
full data deschipti'on of the environment. The simulation requires as 
inpijt a set of surfaces, material conditions of e^ch'^urface, and ttie 
/internal environment of each subset of surfaces'jQ^g a feature. 

^ The general en viro^'"^"* and internal . fe.atur\e ^Environment provide • 
boundary conditions for each surface for several Simulated hours prigr 
to the image ti^e. The complete computation carrfce either executed each 
time surface shades are updated .or stored in tables i ' Either way, th'e^.same 
demands are placed on the scene dats^aje- \^ 

/ the york irj Sectior> IH shows v^^^^rfacp c^^ss section.^. sbsorptfV- 
ity, orientation, and environment all pl3V significant roles in.formfng " 

sigfmture-. These'data are needed in anr potential algorithm that 
generates accurate tdneS- ' 

A high deghge of , ^simulation Accuracy with rfSpect to tonal rendition^ 
does not appear to be hecessary for the training application.. However, 
for passive infh^red sen/ors, the intens^ities must be based upon.param- 
eters which inclydt ^prface cros^ section, absorptiyi'ty, orientation and 
environment. - ' 

/The. number of thermal Paths (conductive, convective, and radiative) 
an3~the time integration of their effects on a surface have a major im- 
jpact on data base requirements. However, all paths, except to the sun, 
can be ignored (extreme simp! ificatron) provided the surface cross sect^n, 
solar absorptivity, sensor band emissivity. orientation, and solar envir- 
onment are known. ■ , ' ^ 1 ^ ' 



Scene Data 

VSheir . ^ 

^ - set of surface? 
Surface Materials 

- solar absiorptivity 

• il^'ouhd/sky absorptivity 

; - sensor band absorptivity 

• surface cross section 
Feature Conditions ^ 

- internal thermal environment 



goviromnent Data 



Sun 

- portion 

- direct radiance 

- diffuse radiance 
Sky 

. - thertnal 
^Ground 

^ thermal radianfce 
Atmosphere 

- air temperature 

• wind vector 

• attenuation rate 

- precipitation 



radiance' 




j^im:^P^-Oata Base 




Simulated Mission 
Conditions 



(From Standard Type 
Environments or Actual 
Measurements) 



Sensor Data 

Sensor Band 
Field of View 
^caiyformat ' 
Point Spread Function 
Response Curve 



Sensor Performance 
Specifications 



Figure 28. Categories of Simulation Data and Their Sources forVlIR 
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5.1.2 LLLTV - ^ . ■ ' . 

. — : — ■ ■ ■..) _ 

" • " TRe si^^^^^ of "passive visible isen^^ There 

■ , are choices as to casting*5hadows and type of reflection (^"ffuse, diffuse 
* ^lus specular, or general' bidirectidhaUreflection) , but these chcnces .' / 
'/'have little impact on data base requirements. Figure 29 presents «ite- 

gories af data -and expected sources. The demands are much less than for 

FLIR sinyjTation. 

t 

The most significant requirements are qn surface material paraoeters. 
Sensor band, reflectivity anri surface condition have important impact on 
surface radiance. In the visible band (^ctually a larger band from 0.2 / 
. to 3 microns), pigments and many other compounds have a- significant and 
high,^ spectrally dependent effect 6\\ reflectivity. 

5.2 AVAILABLE PARAMETERS IN THE^DMAAC DATA BASE ' 

The^MAAC data base serves as a scene data base meaat to be comple- 
mented with environment and sensor data bases suppl iedelsewhere. It was 
basically designed to serve radar simulators. 

^ The DMAAC digital data base provides a ground shell in the form of a 
terrain file and a planimetric (culture) file. The basic resolution^is ten 
or more times lower than a typical E-0 sensor, but the identification of 
discrete features, especially unique significant features, and the use of, 
sta^istiCi|J measures, such as percent tree cover and number' of structures 
per square mile, increase the shape dfifta available to model a shell. 

The manner in which the shell is represented as a set of feature 
gi|ound perimeters produces shape distortions of walls, -roof s,jand small ^ 
structures. The effects of shape distortions and omissions are notice- 
able, as seen in Section IV, but the shapes^are generally i^iterpretable 
in E-0 imagery. ^ ^ 

Surface matisrial paramete^j^must be derived from the material code 
.and feature ID codeK)f each feature. The correlation of surface material 
parameters with the given data^ is loW. For example, roofs and wall's are 
mt distinguished, and surface- condition (smpoth or rough), surface 
reflec^vity in the visible band, surface^. Cross section, etc., are 
nat weln correlated with availa||fie data, "^^^i"^^ ' 

^ • • ■■" /..''■ ■■ . 
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Categories of Data 



Sources 



Scene Data 



Shell , 

- ' set of surfaces . 
Surface Materials 

- - sensor band reflectivity 

- surface cen^it ion (i.e., flat, 

corrugated) 
Feature Conditions 

- internal lighting^ 



DMAAC-Like 



Sce?le Data Base 
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Si>n 

- position 

- ^direct sensor band radiance 

- 'diffuse sensor band radiance 
Atmosphere 

- attenuation nate 



Simulated Mission 
Conditions ^ 



Sensor Data • f 



Sensor Band- 
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Field^of View - ' / . \ Sensor Perf5nnance 

Scan Forniat * 
^ Point Spread Function > { Specifications^ 

Response Curve 




. Some feature conditions"; such as internal environment, are correlated 
with the feature ID jcode so it can be extracted. « 

In summary, much scene data is available but only in' forms correlated 
with the physical parameters, of the scene and 'its features. Some important 
data, such as" visible band reflectivity and other surface parameters, are 
not included in the DMAAC data base. ' 

5.3 RESOLUTION BETWEEN NEEDED AND aVaILABLE DATA 

TheUon.al model developed i*n thks effort reguires more data than are 
presently available in the DMAAC dsJty^ase in order to per/orm acceptable 
tonal assignment. ' The cultural file of the DMAAC source data contains 
information which reflects the general architectural character of the 
rejjfesented features. Although the provided resolution is low compared 
to the capability pf FLIR -and LLLTV sensors, this" file remains thfe best 
source of real-world data available for use in the development of large 
cul tural . data bases. \ 

It is possible that this file could be expanded tp provide more of 
the- data .necessary for sensor simulation^ Such an expansion would, in- 
clude a more detailed description of the geometrical and material prop- 
erties of a feature as, well as visibly b^i^d reflectivity, surface 
cross section characteristics, and other -relevant parameters. 

, There is, of course, .considerable latitude at this, point and continued 
studies should identify idditional deifnands which can be placed on the data 
'b^fise as opposed to those wkich of .necessity must be imposed on simulation 
contractors. - x ' r ' - 
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■ VI. CONCLUSIONS AND RECOMMENDATIONS 

;.. . , r ■ ^ ... 

In this program, a model of the thermal properties of material s^ 
based on theoretical thermo-electromagnetic models was derived. Alsa, 
a data base of the dense, cultural hospital scene was constructed according 
to the DMAAC specifications. In .addition, a program to evaluate the tonal 

modfl and generate imagery of the scene was designed and implemented. 
The simulated imagery shows the complete diurnal and seasonal varia- 
tion of the LLLTV and FLIR sens&r imagery. A key program feature is 
that the LLLTV and FLIR, imagery was generated for the/xhree DMAAC 
levels: lA, HA, andyJilB. Finally, an initial assessment was mac^e 
of the impact on DMAAC to include in the data base the additi'tfna^r des- 
cription of the scene which is necessary to generate the passive *ima^ery.^ 
However, the recommendations of adding the required infonnatipn have been 
of general nature and contain no specific data-^definitioh. Thi^ is 
because the data requirements of the e>ci$ting ^model ar^' rather large. 
Thus, it'-'remains to be determined what^lesser amount of iriformatlon' is 
, required to generate the images for different user applications* . . 

Thus, sensitiv-ity analyses s-f^uld be performed to test the effect 
of the many input parameters ion the resultant- imagery. Then factors of 
Tow-sensitivity can b& made constant lind removed from the model. Factors . 
which have the same sensitivity cambe grouped together. Also, some" in- 
put parameters may have an effect ^YP^he- imagery when they are present 
together and little or no effect when- only one "element is present. , . - ' 
These trend analyses should greatly Vgduce the input c|ata requirements. 

^Once the data requirements have b'een refined, the next step is to 
determine which data should be inc.ludecf iri the data base and how they can 
be included. . Only those factors that affect the DMAAC data b^.e should y^j^ 
be identtfTed, since a number of modeT pgrfamlters are concernedjwith the ; 
descriptioQ of this , sensor or environmenfe'^3^es^ of. this stu^y wiH 
be a specific set of data base descriptors .arV^. the correspCH\ding require- 
ments^ On the DMAAC'data base. 
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j|P^e tonal modeling and its subsequent incorporation withthe geo- 
mAric data base iVito a computationally efficient routine provide a 
practical yet realistic method of computing ttie tonal ^operties of 
cultural scenes. However, the model needs to be augmelrod in order • ) 
to calculate the tonal properties of commonf^ natural features. A study] 
effort for^his purpose would capit^alize on the existing theoretical 
model; thefffnite difference equation solution, and the'software ^ 
implementation. The effort is to model the thermal properties 
of natural features. * . - t - r • 

The work reported here was all directed toward generating images 
0^ the hospitaj scene from the AFHRL Seasonal Sensor Handbook. There 
are other scenes in the^^hanjdbook which have merit from a tactical or 
training context. In fact, LLLTV and*FLIR imagery could be generated 
for any\)f the 1^ scenes contained in the AFHRL Seasonal Sensor 
Handbook.' In addition, scenes from-existing or intentied data bases , 
used: for/'pilot- traifring or scienes of a^tactical interest for bomber/ y 
navigator training could also be used.' Thus, the principles derived ^ 
here could be applied to dther sefiTies to generate their passive 
signatures. 
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